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We study the DC spin current induced into an unbiased quantum spin Hall system through a 
two-point contacts setup with time dependent electron tunneling amplitudes. By means of two 
external gates, it is possible to drive a current with spin-preserving and spin-flipping contributions 
showing peculiar oscillations as a function of pumping frequency, electron-electron interaction and 
temperature. From its interference patterns as a function of the Fabry-Perot and Aharonov-Bohm 
phases, it is possible to extract information about the helical nature of the edge states and the 
intensity of the electron-electron interaction. 



PACS numbers: 73.23.-b, 72.25.Pn, 71.10.Pm 
I. INTRODUCTION 

In recent years the edge states of two dimensional topo- 
logical insulators, showing the quantum spin Hall (QSH) 
effect, have been subject of attention from the condensed 
matter community^ These states have been theoretically 
predicted in graphcne layers with spin-orbit interaction 2 ^ 
and strained semiconductors^ and experimentally ob- 
served in Mercury- Telluride quantum wells 5 -—. Their 
more relevant feature is the helicity^, namely the fact 
that opposite spins propagate in opposite directions along 
the same edge. As far as time reversal symmetry is pre- 
served, the backscattering between these different spin 
channels is indeed suppressed £ 

Various proposals have been done in order to ex- 
tract information about the helical nature of these states 
by means of transport measurements. In the simple 
quantum point contact (QPC) geometryiS signatures of 
helicity and electron-electron interaction can be found 
in the power-low behavior of the spin current as a 
function of the source-drain voltage in a two terminal 
configuration^ The presence of a single peak of the dif- 
ferential conductance in an extended contact geometry is 
also a clear signature of the lack of the spin-charge sepa- 
ration typical of these systems Moreover, by properly 
tune the voltages in a four terminal configuration it is 
possible to generate pure charged and spin currents mak- 
ing these systems suitable candidates for the development 
of spintronicJ^r— 

Extremely interesting are also the interferometric 
properties of these states. Indeed, in a two QPCs ge- 
ometry, the presence of external gates can lead to Fabry- 
Perot and Aharonov-Bohm oscillations without need of 
magnetic fieldsJ^ This fact, together with the possibility 
of switching the current through the gap in the spectrum 
induced by the hybridization of the states on opposite 
edges^i, leads to great advantages in terms of experimen- 
tal feasibility In this geometry, with all the terminals un- 



biased, it is also possible to induce DC spin and charge 
currents by periodically modulating in time electron tun- 
neling amplitudes or external gate voltages. Such kind 
of pumping currents have been analyzed in details for 
the non-interacting case in terms of the scattering ma- 
trix formalism in Ref. [Til . However the investigation of 
the interacting case is still an open problem and requires 
a completely different approach based on the helical Lut- 
tinger liquid mode h 11 ^ 3 

Aim of this paper is to carry out this analysis on the 
spin current. We will consider a two QPCs geometry in 
an unbiased four terminal configuration, focusing on the 
amplitude-modulating tunneling case^, where a DC spin 
current can be induced by means of electron tunneling 
amplitudes periodically modulated in time. We will take 
into account both spin-preserving (SP) and spin-flipping 
(SF) contributions admitted by the time-reversal sym- 
metry. Due to the helical properties of the edge states a 
spin current can be induced only by electrical means, dif- 
ferently from what happens for the spinfull Luttinger liq- 
uids, where a magnetic field is requiredj 19 ' 20 Fabry-Perot 
and Aharonov-Bohm interference paths are induced by 
applying different gate voltages at the opposite edges of 
the sample. We will focus on the first perturbative order 
in tunneling. For high enough pumping frequency the ef- 
fect of the distance between the two QPCs can be taken 
into account in terms of modulating functions showing 
oscillations that strongly affect the behavior of the cur- 
rent. The study of such oscillations as a function of the 
frequency, together with the investigation of the inter- 
ference patterns as a function of the Fabry-Perot and 
Aharonov-Bohm phases, will provide important informa- 
tion about the interaction between the electrons along 
the edges. 

The paper can be divided as follows. In Sec. |H] we 
introduce the helical Luttinger liquid description for the 
edge states of the two dimensional topological insulators. 
We describe the two QPCs geometry in which a DC spin 
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current can be generated by means of periodically driven 
tunneling amplitudes out of phase. For this geometry 
we analyze the SP and SF contributions to the spin cur- 
rent both at zero and finite temperature by taking into 
account the presence of the contacts through modulat- 
ing functions. We discuss the universality of the pumped 
spin charge as a function of the interaction. In Sec. IIIII 
we show the behavior of the modulating functions as a 
function of the pumping frequency for different interac- 
tion strengths and temperatures. We than investigate the 
spin current both by varying the frequency and also as a 
function of the Fabry-Perot and Aharonov-Bohm phases 
induced by the presence of external gates. Sec. |IV] is 
devoted to conclusions. 



II. MODEL 

We investigate a QSH system with a single Kramers 
doublet of helical edge states. Due to their helical 
properties^ one can consider right-moving spin up and 
left-moving spin down electrons on the top edge (T) and 
the opposite on the bottom edge (B) (see Fig [T|) . 

The corresponding Hamiltonians for free electrons 

are 11^3 



can write the electronic operators as 



H 



T(B) = -ifa'F I dx ( : ^R ini) d x ^R, n i) 



L,4.(f)^*i4(t) 
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where ^jj-f (^L,t) annihilates right (left)-moving elec- 
tron with spin up, and analogous for the spin down, and 
vf is the Fermi velocity. With : ... : we indicate the nor- 
mal ordering with respect to the equilibrium state where 
all states below Fermi level are occupied. 

Concerning electron-electron interactions, the only ad- 
mitted terms are the ones that preserve the time-reversal 
symmetry of the system 2 ^, namely the dispersive 
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with <Pr/l,<t( x ) ( a — 1\ I) bosonic fields, -Fr/l^ Klein fac- 
tors necessary to give the proper commutation relations 
between electrons belonging to different edges, a finite 

(T) ( B) 

length cut-off, k F ' and k F Fermi momenta in the top 
and bottom edge respectively. As a consequence of diffcr- 

(T) ( B) 

cnt gate voltages applied to the edges one has k F ^ k F 
(see Fig. [J). 

The bosonic fields <Pr/l,<t( x ) are related to the electron 
density along the edges throug h 22 i 23 
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It is useful to introduce the helical basis for the bosonic 
fields on the upper and lower edge^ 3 - 

<Pt(b){x) = -j= [<Pl,iw(x) - <PR,t(i)(x)] , (9) 

with the canonically conjugated fields 

0t(b)O) = -j= [(p L . w) (x) + ip RMi) (x)] . (10) 

The total low energy Hamiltonian assumes the typical 
form of an helical Luttinger liquid^ - — 
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the interaction strength and 
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locity. In the following we will consider the pure Coulomb 
repulsion (9411 = g-2± > 0) with v = vf/K and K < 1. 



and the forward scattering 



A. Tunneling processes and gates coupling 
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We will neglect possible Umklapp terms, which are im- 
portant only at certain commensurate fillings.— 

Through the standard bosonization techniqu o 22 i 23 one 



As shown in Rcf. 19, it is possible, in a two QPCs 
geometry, to induce a DC current in a grounded one- 
dimensional electron system by means of a quantum 
pump. Tunneling processes are periodically driven in 
time, with a proper relative phase. Here, we consider 
a similar set up for an helical Luttinger liquid in a four 
terminal configuration, where all the external reservoirs 
arc unbiased. 
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Figure 1: (Color online) Top view of the four terminal inter- 
ferometric setup. On the top edge one has right-moving spin 
up electrons (full line) and left-moving spin down electron 
(dashed line). The opposite holds for the bottom edge. The 
amplitudes of both SP and SF electron tunneling at the QPCs 
placed at x = ±o are modulated in time with frequency to and 
phase difference 9. The presence of the top (Vr) and bottom 
(Vb) gates can lead to a mismatch between the Fermi levels 
of the two edges (kp 7^ kp). Shaded grey areas indicate 
the unbiased terminals. 



The more general time dependent electron tunnel- 
ing Hamiltonians that preserve time-reversal symmetry 



top and bottom gates reads: 

H g = dx eVr [-d x (fRf(x) + d x ipn(x)] 

J —a 

+ eV B [~d x tp R i(x) + d x tp L f(x)] . (15) 

As anticipated above, Vr and Vb shift the electronic spec- 
trum and their difference breaks the degeneracy between 
top and bottom boundaries giving kp' B ^ = kp + Kx.B be- 
ing kt.b = KeVr,B /fovF- Note that this definition nat- 
urally reduces to the non- interacting on o 16 i 18 for K = 1, 
as expected. 

The spin current flowing longitudinally in the system 
is defined asi£ 

I spin (t) = ^(N R , t -N L ,{) (16) 
with number operators 

/+00 
dx:*t ta (x)* Vlir (x):. (17) 
-00 

By evaluating the commutation relations of the tunnel- 
ing Hamiltonians in Eqs. (fT12]) - (|13p with respect to the 
above number operators one obtains 




(13) 

being £r l = ±j which represent the SP and SF contribu- 
tion to tunneling respectively. The latter term is strictly 
zero for inversion symmetric systems due to spin conser- 
vation, however gates and deformations of the edges can 
locally modify the intensity of the spin-orbit interaction 
leading to spin flipping processes^ 

In order to generate a DC spin current out of the pump 
one needs time-dependent tunneling amplitudes. They 
can be realized by additional periodically driven local 
electrostatic gates placed at x = ±a 



7j(x, t) = Xj [S(x + a) cos(cj£ + 9) + S(x — a) cos(wt)] 

(14) 

with Xj (j — sp, sf) a real parameter according to the 
requirement of time-reversal invariance, to the pumping 
frequency and 9 a phase factor assumed for simplicity 
equal for both SP and SF term. 

In addition to the tunneling terms, the coupling to the 
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(20) 

the SP and SF contributions to the current respectively. 

It is worth noting that the current in Eq. (|18p was 
already evaluated in Refs. [Tllfl^ with time independent 
spin tunneling contribution only. 

Here we focus only on the spin current, however the 
charged current flowing into the system can be obtained 
in an analogous way by exploiting the relations derived 
in the so called XYZ decompositio n 10 ! 26 that holds for 
helical edge states in a four-terminal setup. 



B. Spin current 

The expectation value of the spin current in Eq. (|T8| is 
given, according to the Kcldysh contour formalism (see 
Refs. [27I - I29I and references therein for detailed discus- 
sions), by 
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where all the operators have to be considered in the in- 
teraction picture with respect to "H e // in Eq. pip . The 
time-ordering along the Keldysh contour C is indicated 
by Tk an d r) = ± labels the upper and the lower branch 
of the contour respectively. 

At the first order in the total tunneling Hamiltonian2S 
one has 



d j (t)) = -- 2 E 
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where the thermal average is taken with respect to "Heff- 
Due to the peculiar time dependence of the tunneling 
amplitudes in Eq. (|15|) . the spin current in Eq. (|23l) is 
characterized by both a DC and an AC component (the 
latter with period 2w).— In the following we will focus 
only on the DC part. It can be written as the sum of SP 
and SF contribution, namely 
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with 
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= 2(4 1) ±fc (2; 



v F j u Fabry-Perot and 

Aharonov-Bohm like phase factorsJ^— 

It is worth to note that a non zero phase difference 
in the time dependence of the tunneling amplitudes of 
the two QPCs (9 ^ 0) is crucial in order to have a DC 
component of the spin current. This is in agreement with 
the general statements of the parametric pumping^ 

The bosonic Green's functions that appear in the Ker- 
nel of the integrals above are22r— 

W R {x,t) = c ( + ) W + (x 1 t) + c ( K ) W-{x,t) (27) 
W L (x,t) = c i k ) W+(x,t) + c ( + ) W-(x,t), (28) 
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are the interaction dependent tunneling coefficients, T(x) 
is the Euler Gamma function^, T the temperature and 
lj c = v/a the energy bandwidth. The latter quantity 
fixes the limit of validity of the helical Luttinger liquid 
picture and can be related to the energy gap between the 
bulk conduction and valence bands of the heterostructure 
which, in realistic experimental setup&i^, is A « 30 meV. 
According to this assumption one has oj c ~ A/S ~ 50 
GHz. 



1. Zero temperature 

At zero temperature, Eq. ([30]) reduces to 

W (0) (i) = -ln(l + iu c t) (32) 



and the integrals in Eqs. (|25[) - (|26|) can be evaluated. The 
SP term can be written as 
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and ujo = vp /2a. 

The modulating function i o 37 i 38 
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with Jp_ 1 (x) the Bessel function of the first kind^, while 
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Pf \E) = 



E 



•9(E) (36) 



is the Fourier transform of the zero temperature elec- 
tronic Green's function vi°'(t) = e^ w<0) W . 

Note that the function also appears in the calcu- 
lation of the spin current in the single contact geometry 
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(both point-like and extended) for a two-terminal config- 
uration with source and drain at different voltage leading 
to a characteristic power-law behavior . n ' 12 

In an analogous way the SF contribution reads 
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where we introduced the new modulating function 

J (0) (£i,6,*) = \ [e-"iF 1 (Ci,a+6;2 i x) 

+e-« 1 F 1 fe,a+6;2ix)] (38) 

with i-Fx(fl) z) the Kummer Hypcrgcomctric confluent 
function^. Note that the above function admits the limit 
J<°) (2,0, a;) = cos a; anc& (£,£,x) = 



2. Finite temperature 

At finite temperature, as long as k B T <C huj c and 
u; c i 3> 1, one has 
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Also in this case it is possible to write the SP and SF 
contributions in a factorized way as 
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The finite temperature modulating functions an 
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; { . . . } indicates the imaginary part and 
2-Fi (a, 6, c; z) is the hypergeometric function^ The finite 
temperature electron Green's function is (for HE, k B T <C 

hhj c ) 
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being $ [a, b] the Eulcr Beta function^ Obviously, ac- 
cording with the above definitions, one has H(^,p,0) = 
H(°Ht,p) and J(d,6,P,0) = J (0) (6,6,P)- 



C. Universality of the pumped spin 



The pumped spin is obtained by integrating the DC 
spin current over a period and one can immediately infer 
that it is proportional to the area in the parameter space 
encircled by the pumping parameters, i.e. -| A 2 sin (9, 
where i = sf, sp. Let us note that in the limit of low 
pumping frequencies lo <C wq (ujq = up /2a), the trans- 
port problem is equivalent to consider a scattering off a 
single impurity localized at x — 0. By using the Kane 
and Fisher argument^ when considering repulsive in- 
teraction K < 1, the impurity is a relevant perturbation 
in the infrared limit and thus in the presence of a small 
bias all the current is backscattered or spin-flipped. For 
the spin pump this picture implies the maximal pump- 
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ing response in the adiabatic limit, i.e. the charge or the 
spin pumped per cycle is quantized. In the case ofspin- 
less particles this has been exactly proved in Ref. [20 for 
K = 1/2 where the problem can be mapped to the time 
dependent scattering involving free chiral fermions and 
an impurity state. 



III. RESULTS 
A. Modulating functions 

Before entering into the details of the behavior of the 
spin current, it is useful to analyze the peculiar forms of 
the modulating functions both at zero and finite temper- 
ature. Fig. [5] (a) shows for different values of the 
interaction strength if as a function of uj/ujq . As already 
discussed in Ref. H3 in the framework of the interfero- 
metric properties of the fractional quantum Hall effect, 
it shows an oscillating and rapidly decaying behavior as 
a function of the pumping frequency due to the presence 
of the Bessel function. These oscillations are related to 
the finite distance between the QPCs. For very close 
point contacts this behavior disappears. The envelope of 
the curves decreases algebraically as (w/wo) _dR (see Eq. 
(|34[) . The positions of the zeroes move to higher frequen- 
cies by decreasing K and consequently the number of 
observed oscillations in the considered range of frequen- 
cies decrease by increasing the interaction. An analogous 
scenario can be found for the modulating function J*- -* 
(see Fig. H](b)) that is related to the SF tunneling pro- 
cesses. However, in this case, the oscillations are only 
slightly suppressed at high frequency. 

The comparison between the zero temperature and the 
finite temperature cases is shown in Fig. [3] for K = 0.6 
(analogous results can be found for other values of the 
interaction). It is possible to note that, for T <C To = 
HuJo/kB, oscillations of H in Eq. (pl2"]) (see Fig. [3] (a)) are 
progressively suppressed by increasing the temperature. 
However, the periodicity and the zeroes positions are not 
affected. In the same temperature regime the function J 
in Eq. P5|) is only slightly affected by the temperature 
with respect to H, as can be seen in Fig. [3] (b), so spin- 
flipping tunneling terms are quite robust. For T ^> T = 
fuvo/ks the features of both modulating functions are 
washed out and additional tunneling processes need to 
be taken into account^ 

Note that assuming a fa 1 /im and vp w 5 x 10 5 m/s, 
as appears from experiment o 6 ' 11 ^ 3 , one obtains wo ~ 250 
MHa^i and consequently To w 2 K. Being the typi- 
cal temperature at which experiments are carried oul£i£ 
T « 30 mK it is possible to conclude that thermal effects 
can be reasonably neglected in the relevant experimental 
cases. This allows us to focus in the following on the zero 
temperature case. 

We also note that at low frequency (oj <C ujq), the effect 
of the modulating functions is unimportant (H, 
Moreover for high enough frequency the SP contribution 
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Figure 2: (color online) Modulating functions (a) 

H (0) (d K , Klu/ljo) and (b) J (0) (2c£° , 2c { ~ } , Ku/u)o) as a func- 
tion of oj/luo for K — 1 (non-interacting case, full black curve), 
K = 0.8 (dashed green curve) and K — 0.6 (short dashed blue 
curve) . 



disappears and the physics is dominated by SF tunneling 
processes. This fact is a consequence of the differences in 
the interference paths induced by the helical properties 
of the edge states 



B. Current 

The evolution of the zero temperature SP and SF 
current contributions as a function of the pumping fre- 
quency uj is shown in Fig. H for different interactions 
and for a peculiar combination of the Fabry-Perot and 
Aharonov-Bohm phases (0fp = <Pab — ?r/4 (mod 2-7r)). 
It is possible to note the emergence of the power-law be- 
havior typical of the helical Luttinger liquid P^d ( w ) oc 
{uj /uj c ) k+ t<~ 1 , already predicted for the current as a 
function of the source-drain bia a 11 ! 12 , but strongly mod- 
ulated here by the presence of the oscillating functions 
if(0) (Fig. 11(a)) and J (0 ) (Fig. 11(b)). 

Fig. [5] shows the total DC pumped spin current. Ac- 
cording to the previous considerations, at high frequency, 
it is driven only by the SF contribution due to the alge- 
braic damping of the SP one, therefore a detailed study 
of this regime can provide important information about 
the role played by this kind of processes in the dynamics 
of the system. 

Even more interesting is the dependence of the cur- 
rent on the Fabry-Perot and Aharonov-Bohm phase as a 
function of the pumping frequency and the interaction. 
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Figure 3: (Color online) Modulating functions as a 
function of uj/ujo for K = 0.6 at various tempera- 
tures: k B T = (full black curve), k B T = 0.2hu 
(dashed green curve) and k B T = hcoo (short dashed 
blue curve). (a) H(dfc, Kcu/ujn, Kk B T/hu)n) and (b) 
J(2c ( + ) , 2c { ~\ Klj/cuo, Kk B T/hu ). 



As already shown in the non-interacting case this leads 
to very intriguing patterns.— This is due to the different 
modulation in amplitude and sign of the SP and SF con- 
tributions to the spin current. Some examples of these 
possible patterns are shown in Fig. [5] for various fre- 
quencies at fixed interaction {K = 0.8). At low pumping 
frequency (top left panel) both the SP and SF contribu- 
tions affect the current leading to clear periodicities both 
in (fipp and <^ab ^ The amplitude of the oscillations cru- 
cially depends on the ratio A s //A sp . Typically it is rea- 
sonable to assume that these tunneling amplitudes have 
the same order of magnitude with A sp > A s /. 16 i 18 

In this kind of setup it is also possible to tune the 
pumping frequency in such a way to have a SF contri- 
bution to the spin current equal zero, while the SP is 
present (top right panel); this leads to a clear period- 
icity only driven by </>fp- The observation of this ver- 
tically striped pattern can give a precise indication of 
the positions of the zeroes of the SF contribution and 
consequently of the value of the interaction strength. In 
the opposite case of null SP contribution and SF differ- 
ent from zero horizontally striped patterns appear driven 
only by </>ab (not shown). This configuration is typical 
also of the high frequency regime (bottom panels) where 
the SP contribution becomes negligible due to its alge- 
braic damping. The differences in the patterns between 
the left and right bottom panels depend on the sign of the 
modulating function J. The strong dependence of the in- 
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Figure 4: (Color online) (a) SP and (b) SF contributions to 

A 2 h ~ 

the DC spin current (in units of In = , .f F — sin 8 and In 

A 2 ft 

A i* sin 9 respectively) as a function of u (in GHz) for: K = 
1 (non- interacting case, full black curves), K = 0.8 (dashed 
green curves) and K = 0.6 (short dashed blue curves). Other 
parameters are: vf = 5x 10° m/s, a = 1 fim (wq = 250 MHz), 
a — 10~ 8 m, 0fp = <Pab = tt/4 (mod 2tt). 

terference path on interaction can be seen in Fig. [7] that 
compares the interacting (left panel) and non-interacting 
(right panel) case. 

It is worth to note that the patterns of interference 
described above can be observed only if both SP and SF 
tunneling processes occurring at the QPCs are present. 



IV. CONCLUSIONS 

We proposed a two QPCs experimental setup involv- 
ing helical edge states of a two dimensional topological 
insulator. Here, a DC pumped spin current is induced by 
periodically modulating in time SP and SF electron tun- 
neling amplitudes out of phase. We analyzed the different 
modulating functions associated to these processes both 
at zero and finite temperature. We discussed also on the 
universality of the spin pumped in a period. We found 
that the power-law behavior typical of the helical Lut- 
tinger liquid expected for the spin current in a single QPC 
is strongly affected by interference effects and interaction 
strenght. We investigated the interference patterns cre- 
ated by Fabry-Perot and Aharonov-Bohm phases induced 
by external gates connected to the edges. The measure- 
ment of their features, and in particular the locations 
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-0.1 

1 2 3 4 5 

u {GHz) 

Figure 5: (Color online) Spin current I sp in(u) (in units of 

A 2 ft 

Io = , -/ p — sin 6) as a function of ui (in GHz) for: K = 1 
(non-interacting case, full black curves), K = 0.8 (dashed 
green curves) and K = 0.6 (short dashed blue curves). Other 
parameters are: vf ~ 5x 10 5 m/s, a = 1 /im (ojo = 250 MHz), 
a = 10~ 8 m, X^f/Xsp = 0.5, FP = 0ab = tt/4 (mod 2vr). 




7T 27T 7T 2tt 

0FP 0FP 



Figure 6: (Color online) Density plot of the spin cur- 
rent Ispin{u) (normalized with respect to I p0 int{ui) = 

A 2 ft ~ 

4n s f a 'i sm6V2d K {ijj)) as a function of the Fabry-Perot 4>fp (x 
axis) and Aharonov-Bohm cpAB (y axis) phases for different 
pumping frequency: u: = 100 MHz (top left panel), uj = 500 
MHz (top right panel), ui — 1 GHz (botton left panel), oj = 2 
GHz (bottom right panel). Other parameters are: a — 1 
/im, v F = 5 x 10 5 m/s, (w = 250 MHz), a = 10~ 8 m, 
X 2 sf /X% = 0.5, if = 0.8. 




7T , 27T JT , 2* 

0FP <PFP 



Figure 7: (Color online) Density plot of the spin cur- 
rent Ispin(oj) (normalized with respect to I p0 int(u) = 

A 2 ft 

47V Y a i sin6V2d K (<jj)) as a function of the Fabry-Perot (/>fp (x 
axis) and Aharonov-Bohm </>ab (y axis) phases for interaction 
K — 0.8 (left panel) and K = 1 (right panel, non-interacting 
case). Other parameters are: a = 1 ^m, vf = 5 x 10 J m/s, 
(cjo = 250 MHz), a = 10~ 8 m, X 2 sf /X% = 0.5, w = 500 MHz. 
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of their zeros, represents an important tool in order to 
extract information about the strength of the electron- 
electron interaction in the system. 
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